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Abstract 

Einstein's equations for a Robertson- Walker fluid source endowed with 
rotation are presented upto and including quadratic terms in angular ve- 
locity parameter. A family of analytic solutions are obtained for the case 
in which the source angular velocity is purely time-dependent. A subclass 
of solutions is presented which merge smoothly to homogeneous rotat- 
ing and non-rotating central sources. The particular solution for dust 
endowed with rotation is presented. In all cases explicit expressions, de- 
pending sinusoidally on polar angle, are given for the density and internal 
supporting pressure of the rotating source. In addition to the non-zero 
axial velocity of the fluid particles it is shown that there is also a radial 
component of velocity which vanishes only at the poles. The velocity 
four-vector has a zero component between poles. 



1 Introduction 

Perturbation techniques, so important in the General Theory of Relativity, have 
frequently been applied successfully in the description of slowly rotating compact 
perfect fluid sources. Of considerable importance has been the analysis of Hartle 
P| who presented the equations for the equilibrium configurations of cold stars 
up to and including the first order of angular velocity parameter. This work 
has formed the basis of an extended first order analysis for example, Kojima 2 
and [3| Beyer & Kokkotas 0], Abramowicz et al [5], to the address the problem 
of the non-radial quasi-periodic oscillations of rotating compact sources and the 
resulting r-mode spectrum of relativistic stars . 
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In examples of non-compact rotating sources Kegeles [H] and Wiltshire [7] 
successfully applied the perturbation method to a slowly rotating non-equilibrium 
configuration. In this case a Robertson- Walker dust source endowed with rota- 
tion was successfully matched to the Kerr exterior solution of Einstein's equa- 
tions to the first order in angular velocity parameter. 

Recent second order perturbation analyses have largely been confined to non- 
rotating cases. For example, Salopek et al [S], Russ et al [5] use the method to 
discuss the gravitational instabilities of an expanding inhomogeneous universe. 
However, there is seemingly an absence of literature on the use of second order 
techniques to describe rotating bodies. This is perhaps surprising, since such 
analyses can reveal the relativistic effects of the spatial distribution of a rotating 
fluid as characterised by internal density and pressure. Such effects are not 
revealed in first order approximations where fluid density and pressure are shown 
to be same as for a non-rotating source. Moreover the second order effects 
result in calculated deviations of spherical symmetry in the fluid boundary of 
a compact or extended body which may be used in the context of the matching 
problem with a Kerr vacuum. 

It is the aim here to present an example of a second order perturbation 
analysis applied to an extended rotating source. In particular a Robertson- 
Walker source will be endowed with rotation. As will be seen the approach 
naturally brings about solutions of Einstein's equations which exhibit the non- 
homogeneities in internal density and supporting pressure due to rotation. 

In the following a non-rotating source will be described using the Robertson- 
Walker metric in the form: 

d<J 2 RW = d v 2 - R 2 (ry) (^^2 + + ? ( 9 ) W 

where R = R(rj) , k = — 1, 0, 1 and where the speed of light c = 1, and gravita- 
tional constant G = 1. The homogeneous density and supporting pressure will 
be denoted by and prw respectively. The fluid source will be endowed 

with rotation which will be characterised in terms of an angular speed parame- 
ter, denoted by q, and the mathematical analysis to follow will be accurate up 
to and including second order terms in q. 

The rotating fluid source will be described in terms of an extended form of 
the Robertson- Walker metric which is taken to be: 

da 2 = (l + Qg 2 ) drj 2 - ^ ^\ (l + lV) df - 2J^R 2 q 2 d^dr 1 

-fR 2 (l + Vq 2 } d0 2 - £ 2 R 2 sin 2 (0) (l + Wq 2 ^ d(/) 2 

-2£ 2 sin 2 (6) R 2 q {Yd£d(f> + Xd^dij) + O (q 3 ) (2) 

in which each of the functions, J, Q, U, V, W depend on £, 9 and 77 whilst X, 
Y depend on £ and r\ alone. Note the this form of metric is an extension of 
the linearly perturbed form of the Robertson- Walker metric used by Kegeles 
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6 . The components of the fundamental tensor from are determined using 
Einstein's equations for a perfect fluid written here in the form: 

G a b =-8nT b a , T b a = (p + p)u a u b -S a bP , (3) 

where p, p are the respective rotating source density and supporting internal 
pressure and u a are the components of the velocity four- vector with the property 
that u a u a = 1. 

The use of seven as yet unknown functions in J5| will naturally gives rise 
to ambiguity in the solution of Einstein's equations. This ambiguity can be 
removed by further specification of the gauge in which the solutions are to 
be determined. Although convenient gauge choices are much discussed in the 
literature, including recently, Bruni et al |l(Jj it is convenient at this stage to 
continue the solution process without further specification of the gauge other 
than that which is explicit from the metric choice 

For brevity the term O (g 3 ) , included in W1 U be omitted from all future 
expressions. The fact that all terms including q n , n > 3 are taken as negligibly 
small will of course be implied. 

In summary it is the aim in the following to determine solutions of Einstein's 
equations in the form J3J for rotating sources described by the metric © 

2 Solution Approach 

Direct calculation of the components of the Einstein tensor for (J2J show that the 
components G\ and G\ are identically zero and so for a rotating perfect fluid it 
follows that the velocity four vector component u 2 = and one form component 
u 2 = 0. It thus follows that the conditions T? = = and T| = = T 2 4 must 
hold upto and including terms in q 2 . Moreover Tf + p — 0. 

In addition the perfect fluid © must also satisfy the following consistency 
relationships 

(Tl+p)(Tl+p)~TlTf = 
(T}+p) (T£+ p) -TlT* = 

(T!+p){Tt+p)-TlT± = (4) 

However by direct calculation G\, G\ , G\ and G| depend on linear terms in q, 
whilst G\, G\, G\ , G\, G\ and G\ depend on quadratic terms in q . Also when 
q = each of these components is zero. It follows that for the Robertson- Walker 
source endowed with rotation that solutions of Einstein's equations must satisfy 
the following perturbation equations: 

Tl+p = 

T 2 +p = 

(Ti+p)(T*+p)-TiTi = 

T3 1 = If + 
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Tf = = T\ 

Tl = = T 2 4 (5) 

The first or second of these equations may be used to calculate the internal 
pressure p whilst the density p is calculated using 

p = T« + 3 P (6) 

where the repeated index indicates summation. It is perhaps worthy of note that 
although the velocity four vector u RW — for the standard Robertson- Walker 
case, the system of equations © do not imply that this condition is retained 
for the rotating source. 

The angular velocity of the source will be denoted by L (£, 77) where: 

Since the fourth of conditions © T3 = may be solved immediately to give: 

Y„ =X, + , (8) 

where h(£) is an arbitrary function of £, it follows from (JJJ) that the angular- 
velocity of the source is given by: 

gvg ~ qX (9) 

16< R 5 (p RW +prw) 

where the suffix denotes a partial derivative. The density Prw and pressure 
Prw for the standard Robertson- Walker metric Q are such that: 

8^ [pnw + P RW ) - + 2 -§ + % (10) 

Moreover, a particle moving in the field of @ will have zero angular mo- 
mentum whenever U3 = 0, so that the quantity: 



U3 = qsm 2 9y/l - kfh^ 
u A IQ-k^R? (p RW + prw) 

will also be zero for such a particle. It follows that the induced angular velocity 
fi (£, ??) of the inertial frame is given by: 

n f (Z,ri) = -qX , (12) 

and that the angular velocity of a particle moving in the field of is: 

n p (S,r,) = p/lEKh (13) 

16tt£ R 5 (p RW +prw) 
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Clearly therefore it is the nature of h (£) which determines the actual angular 
velocity of the system and that h = defines a non-rotating source and that 
the solution of © is then: 

X = $„ Y = $ c (14) 

for some $ = $ (£, 77) . 



3 Simplification of the Perturbation equations 

Using the second equation p — — Tf the first and third of (J5J with © become: 

G}-G 2 =0 (15) 

(G 3 3 - G 2 ) (Gj - G\) - G\G\ = (16) 

In addition the fifth and sixth of © namely, Tf = = and T 2 = = T 2 4 
will be satisfied by taking: 

G\ = (17) 

Gl = (18) 

These equations may be simplified somewhat by firstly, defining a new time 
variable r (rf) through: 

-W = /f (19) 

and, secondly by introducing the functions S (r) and T (r) expressed in terms of 
the density Prw and pressure prw for the standard Robertson- Walker metric 
as follows: 

n I \ 1 / \ 2R TT SR 2 2k rt\r\\ 

^ T > = wf = ^ PRW + Prw ' = — W + ~W + 1¥ ' 

Thirdly the equations (|15|l to (|17(l may be rendered independent of W with the 
aid of the following substitutions: 

J = XY sin 2 9 



U = (1 - kS, 2 ) Y 2 sin 2 9 + U + W 

v = v + w 
w = w 

Q = -fR 2 X 2 sm 2 9 + Q-W (21) 
where each of U, V, W, Q and J are again functions of £, 9 and r. 
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It follows from the transformations 121|) and l|19fl that the metric @ may 
now be written in the form: 



? = ( x + ( Q - Y ] g2) * *> - (1 + f + w lf ] <e 



2Jq2 f dtdT-(l + (V + W)q 2 )d8 2 



-Wq 2 sin 2 (0) dcj) 2 - sin 2 (9) (d(j> + qYd£ + qR 3 Xdr) 2 } (22) 

Note that the entity d(f> + qYdi; + qR 3 Xdr is itself an exact differential only in 
the non-rotating case when equation 114|) applies. 
In this way the first of equations (|f 5J) becomes: 

-£(l-fc£ 2 ) ( Us + V c )- e{UrT R ; VTT) +2(U-V) 



, (V e + Ue)cose h 2 sin 2 9 _ n 
In addition the second of equations l|16[) has the form: 

e (1 - ki') V C£ - 3kfV e + 2(V e - S» + Ult + Q„ 

_<P. + Q.)°°-'_^,(W-f)T + »_) =B (24) 
sin0 ^ 2£ 4 C 4 i? 4 J 

The condition (|f 7|l that G 2 — becomes: 

Only the remaining equation 118fl explicitly contains 14^ as follows: 

2iiV (W, - Q,) + + _ ^ + _ Kcosfl = Q 

it v smO 

Equations (|23|) to I26[) are the hnal forms of the perturbation equations which 
determine U, V, W, Q and J for given h (£) for the metric (|22|l . 

The internal supporting pressure, calculated using the second of J5| is: 



6 



(W T -Q T )R T k (U + 2Q - W) _ (2W TT + U TT ) 

R 7 + R 2 2i? 6 

(1-fcg) ( 2Q i -U i \ Q e \ 

2R 2 \ £ + ^ ? V 2£i? 2 J 

(j 2 (Q 9 + t7e)cosfl _ /iVsin 2 
2£ 2 i? 2 sin 9 4£ 6 i? 6 

whilst the internal density, calculated using 10 is 



+ q ^ru^ _. 1 ^ L + &KpHw (2?) 



8irp = q 




(W ee + Uee + Qee) 

e 2 i? 2 



i? 6 i? 7 v y R 8 R 2 



£ 2 i? 2 ' i? 2 V"'^" 5 £ ' ^ S '' £ 2^ 
v/l - ( 2^R T 6JR T 



g 2 {Qe + Ue + 2V e - 2W e ) cos 9 h 2 g 2 sin 2 9 
2£ 2 i? 2 sine 4£ 6 i? 6 



8ttp w (28) 



Further note that the application of the transformations (|21|l and (|19|) has 
enabled each of the equations (|23|l to l|28|l to be written in a form which is inde- 
pendent of both X (£, r) and Y (£, r) but explicitly contains terms in /i (£) .This 
is expected from <|12ll since one would not expect internal pressure and density 
to be dependent on the frame dragging term X but rather on that function 
defined in © and which directly determines the true angular velocity of the 
source, namely h. It follows that further analysis may continue without further 
detailed specification of X and Y, only h needs consideration. 

Finally, the two equations (|23l) to H25fl contain the four unknowns J,Q, U, 
V and may be used to determine families of rotating extended sources for a 
range gauges and or physical conditions. As a particular example of a solution 
procedure it should be noted that in cases when V is known explicitly then 
equation $Z%$ may be integrated directly to determine U + Q. Thus if either 
U or Q is known then the equation (|25|) may be integrated immediately to 
determine J. 

Finally, the velocity four- vector component u 1 ju A may be calculated through 

U - U (29) 



u 4 T£+p 

and so up to and including quadratic terms in q 2 it follows that 



u 



7 4 



l4 8f(p w +Prw) 



(30) 
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where G\ is given by: 

R5 °l Urn , 1 fVr rr\ , R r 



V^~W IV 2 / € V 2 7 ^ 



^ + + 7 hr - u r )+ -w (Wc -QO Wi-K 



,2 



* J '2^J+^| (31) 



2£ i? 4 T s 2£ sin 6 
Note also that in a similar way ui/u^ may be found using; 

ui = Tf = G\ 

m T^+p 87r (p RW + prw) 

where G\ is 



(32) 



„ 4 _ /i g g 2 sni 2 0Yy/l-k( 2 2 f 2W r€ + Kg | 2Ef T - K 
G i _ ., ,-2 1 



2f i?3 y \ 2i? 3 2^ i? 3 

R T (Qs-W s ) 2Jy/l-k? / J 9e + 4 J \ | 



R 4 tR 3 \2Z^/T~~k?R\ 

q 2 cos 9 Jg 



- (33) 

2 sm6£y/l-k^ R 3 

Thus only l|33[) depends explicitly on Y (£, t) and is a further manifestation of 
the frame dragging effect expressed in (JSJ. This may be removed by choosing 
Y = {). 



4 Characterisation of angular velocity leading to 
analytic solutions 

Consider first equations the three equations 1|23|) to 11251) and note that particular 
solutions may be found in principle by setting: 

U{W,t) = ^sm 2 9 + u 2 (34) 
V{iAr) = ^p-+u A (35) 



J(tO,T) = ^l-ke(u 5 sm 2 d + u 6 ) (36) 

Q(£,6,t) = u 7 sm 2 9 + u 8 (37) 

where Ui are functions of £ and r alone. It is straight forward to show that 
equation (|24|l may be used to determine u-j and U4 since: 

s 7 2i? 4 2 ^ J « 2 y 

9fcf 3 li3, „ , v 

+ - * - I2kfu 3 - 5^3, + 15u 3 - ui (38) 
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and U4 = 0. Moreover direct substitution of these relationships into H25[) shows 
that this equation is satisfied provided that U5 is such that: 

-p- = *ii? 4 T + ^ 2 R A + * 3 (39) 

OT 

where ^1,^2 an d ^3 ar e fully defined in terms of h, u\ and U3 and are presented 
in the appendix. The equation lp!H| and (|3*9f together with w 4 = may then be 
substituted directly into equation (|23[1 to produce a very lengthy relationship, 
reproduced in the appendix, having the general form: 

* 4 sin 2 6 + * 5 = (40) 

where "J 4 and ^5 are again fully defined in terms of h, u\ and 113 and their 
partial derivatives with respect to £ and r. 

Inspection shows that the equation ^ 4 = may in general terms only be 
solved numerically for Mi(say) in terms any given U3(£,t) and /i(£). Remark- 
ably, however there is at least one case where considerable simplification is 
possible, namely when: 

h(o= n e (4i) 

with constant n, for which a basic separation of variables approach reveals that: 

ui (£, r) = n 2 e° (V-i + M 2 ) "3 (£, t) = n 2 ^V 3 (42) 

where tp 1 , ip 2 an d ^3 ar e functions of r alone which satisfy: 

lOJfe (1 + ib» ) + 1pn 
125fc 2 T * I^d Iliz. _ 24fc + lOA;^) = (43) 

- ^ + 7 ^"^ +5 + 4. (49^ - 2 A ) + 14^ 2 = (44) 

Although these equations need to be solved numerically when k = 1, — 1 the 
case A: = yields the analytic result that: 

7^ 3tt " ^i TT +5-0 ^2-0 (45) 

In the following the family of solutions presented will be based upon equation 
Q45J1. Thus using i|41[l and (|19|l substituted in (|13|l and the angular velocity 
flp (£, t) of a particle in the field of (|22|l is: 

Vp(Z,T) = - - Hq r (46) 

16ttR 2 (p RW +PRw) 

It follows that the angular velocity is purely time dependent. The frame drag- 
ging effect (fT2l is 

n f = ~qR 3 X (47) 
where from JHJ and (|19(l X (£, r) is determined through: 

i? 3 X ? = Y T - < (48) 



9 



In cases when Y = so that the frame dragging effect due to rotation is removed 
from <|33l) then 146fl becomes: 
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n, = =f (49) 

where for choice X (0, r) = 0. 

5 Development of solutions with k=0 and with 
purely time dependent angular velocity 

Thus using the solution procedure outlined in the previous section with (|45|l it 
is straight forward to show that the corresponding solutions of (|23|) to l|2tj[l are: 

U(£,9,t) = n 2 V>i£ 4 sin 2 (9 + n 2 Z2 (50) 
V(U,r) = t(^-^-)t^ 2 (51) 



7 V 2 

J (£, 9, t) = n 2 z 5 sin 2 6 + n 2 z 6 (52) 



2 , 4 /3^ 25T 25r 2 \ ? (V Vt + 2) . . 2 2 _ 



Q ({, <?, r) = n 2 " -)- V ;^ fl4 ' | sin 2 + n 2 z 8 

(53) 

^r) = ^^ + ^ (54) 

where z$ (£, r), Z6 (£, r) and zg (£, r) are determined in terms of ip 1 (r), z 2 (£, T ) 
and (£, T ) through: 



9 2 z 6 



d£dr H 4 14 7 J 2£ 2£ 2 £ 



9zg 

97 



' 28 



45t 2 


5^i \ z$ u 


14 


7 J 2£ 


z-? 
4- — 




2£ 




25T 


25r 2 \ £ 6 i? T 




7 J 14i? 4 




R 



(57) 
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The function Ziq (£, t) is arbitrary. The supporting internal pressure (|27[) is 
given by: 



8np = 8 7 rp w +nVsin 2 0|^ + Ii |^(25r 2 -3V. 1 ) + 2^ 
+ 2 8i? io r l^i TT + ^+ 8i?6 

Zg Zg 3z S R 2 Zg T R T ZlQ T R T 



-n 2 q 2 



i? 6 T 2i? 7 T i? 8 i? 7 
£R 2 £ 2 R 2 R 6 R7 



(58) 



with: 



1 3R 2 , , 

8ttpw = _ - — g- (59) 



i? 6 T i? 8 
and the internal density is (|28|l is 

{ ( 25T \ £ 4 i? 

8np = 8np RW +n 2 q 2 S m 2 e^i [ -— + 5T 2 +6iJ 1 j-^(5T + yj lT ) 



Rt u n n 1 / 2z 9? 6z 9 \ £ 



2i?7 VS~o s , ~~oj > R3 y £ ' £2 z 9« y 4i?6 



,4 



i? 8 1 i? 7 V s 4 1 i? 



with: 



3i? 2 

Stt^w = -j^P (61) 



Also using and (|31|l it may be shown that: 

1 u 1 2 2 { &6 . z w Ti - zs s R T 2z 5 - z 2t 



RT u 4 " ? | R 4 T R ' £ 

+n q sin | - (z 5 « - ^ J + *« - ^ ^ (J2) 

Inspection of (|62() shows that in general terms u 1 (£, 0, r) although it is possible 
to choose zio so that u 1 (£, 0, r) = 0. Thus u 1 = is possible for a particle 
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moving in the field of (|22H only for an observer situated on the axis of rotation. 
In this case: 

6 A subclass of solutions 

Consider now the particular case: 

ipl = -t 2 z 2 = UJ2S, 2 z 8 = w 8 £ 2 (64) 
where tt>2 and ujg, are functions of t. The equations l|55|) to l|57|l become: 

^ 5 = ^ 2 ^^ 4 ( Z f + 10 ^) < 65 > 

/ f 2 ^ p4 /25T 5r 2 \ U2rT 

Moreover 

z 9 = a 2 e + a 4 ? 4 (67) 
da 2 3ip 5 R da 4 3ri? / 25T 2 \ 

— + 4t J (68) 

where s/> 5 , "0 6 , a 2 and «4 are functions of t alone. The condition (|63J) is: 

«io = /? 2 ^ 2 + P^ 4 (69) 

^2 n , n (u 2t , \ d(3 A i(j 6 

dr 

where (3 2 — /3 2 (r) and /3 4 = /3 i (r) 



ad /3 4 = /3 4 (r). 
The supporting internal pressure is given by: 



Sirp = 8irp RW + n 2 q 2 sin 2 d{e(^ + ^ 



\8R 2 2R 7 T J 

5 \8R 6 R 6 T 2R 7 T) 
2 2 2 / 75T W8 (3 2 3u> s R* T R T ( /3 2t 
~ U q ^ Tii? 2 " _ Wf + 2tfT + + W P T + IT 

10r 2 w 2rr /? 2t 



i? 2 2i? 6 R 7 

{ 04 Pa. 

\ 2R 7 T R s R 7 J ' R 2 



" n <? f ttt^t; + —5^ - "57 f + "5T { 2ws " ^ ( 71 ) 
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whilst the density is: 



- {2w 8 - 


^2} 


f3o;2 


4a 2 


1 i? 2 





4i? 6 V 

2 2 2 f ZujzRI (uj 2t -10^)R t , 3r 2 4a 4 

2 2 J 3w 2 4a 2 \ , _ s 

+Hq (72) 

Notice that in region of £ = that both the expressions for pressure and 
density are well behaved and in particular when £ = the pressure and density 
of the centre of the rotating source are 

8irp (0, r) = 8np RW + -^§- {2to a - u 2 } (73) 

8np (0, r) = 8np RW + +n 2 q 2 j ^ - ^ J (74) 

Note that the centre of the source will be non-rotating 8np (0, r) = 87ip.Rw and 
87173 (0, r) = 87179^,^ when: 

" 2 = 3i? " 8 = 3i? (?5) 

For this solution (|62|l defining it 1 becomes: 

^ = „ V {W S (|-A) (T6) 

7 Robertson- Walker dust endowed with rota- 
tion 

As a specific example of the above subclass consider the case of dust endowed 
with rotation so that: 

R(r) = *t§ = ^ ( 77 ) 

where k is constant and suppose also that 87175 (0, t) = 8ttprw and 87179 (0, r) = 
8?r Pw in GSJ and Oj tnen to (SIJ| become: 

C/({,0. T ) = - 207 ^" V - n V^si,»0 (78) 
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272 



5n 2 g 2 sin 2 <9 I 159t*£ 2 29£ 4 I 7715r3n 2 g 2 £ 2 

8np= i6. 2 hr + w + 272. 2 (83) 



3r 6 K 6 4k 2 I 136 T 4 K 4j 272k : 

and 

w 1 _ 645k 6 t 7 ?iV£ 3 sin 2 i 
~~ 272 



n . 139r 2 n 2 C 4 . , „ 207 K 4 rfn 2 £ 2 

Q{e,e,r) = — sin 2 6» — (81) 

^ vs " ' ' 16 1088 v ' 

„ r „ A , n 2 K 4 r 6 sin 2 /945r§^ 2 121^ 4 \ 837 K 4 rfn 2 £ 2 

W *' T) = 64 - ^ J K^T^ (82) 

and (|71|l . (|72|) and l|76|l become respectively: 



« 4 , 3nV , 2o U5rH 2 £ 4 2915rln 2 g 2 f , 
8?r /> = XXIfi + XT sm 8 S —^T* XI f + T^TT-o (84) 



(85) 



8 Conclusion 

Einstein's equations for a Robertson- Walker source endowed with rotation up 
to and including quadratic terms in angular velocity parameter have been pre- 
sented. It has been shown that a family of analytic solutions of the equations are 
possible for the case when k = and the angular velocity of the fluid is purely 
time dependent. The corresponding density and supporting internal density are 
explicitly presented in a form containing perturbations from their respective 
Robertson- Walker counterparts. A subclass of the solutions merges seamlessly 
with the Robertson- Walker source at the origin. The work presented here is very 
much a preliminary investigation and further research is now being conducted 
since it is possible that further mathematical analysis will reveal new analytic 
solutions of Einstein's equations with for example spatially varying fluid angular 
velocity. Moreover a numerical analysis will be applied to reveal a broader range 
of properties of the perturbation equations for varying forms of h = h (£) and 
for k = 1,0,-1. Clearly, the determination of rotating sources incorporating 
graviational radiation would also be an important development. Further work 
will also address gauge issues since the analytic solutions presented in this pa- 
per have the property that the velocity four vector components have the general 
property that u 2 = and u 1 ^ except at the poles where u 1 = 0. It would be 
interesting to consider solutions for which u 2 = and u 1 = for all values of £, 

9 and r. 

Whilst perturbation analyses are of considerable importance in the General 
Theory of Relativity there is no doubt that the major goal for future research 
must be in the determination of exact solutions of Einstein's equations for rotat- 
ing sources with physically realistic properties. According to Bradley et al 1 1 1 j 
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there is currently an 'embarrrassing hiatus' in the availability of such solutions. 
It is the hope that the perturbation analysis presented above may provide a 
signpost which lead to the possible discovery of such solutions assuming that 
they exist. 
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A Expressions terms used in section 4 



In the following: 



= ^Ji-te (ao) 



so that: 



h f h ff P 2 3 hi P 2 h 2 

*i(€,t) = * K ^ ^ (Al) 

v ' 2£ 7 4£ 8 2£ 8 

T 2 — 7 T a Q T 



2£ 7 e 8 £ 9 £ 10 

| 3u 3(( 29u 3( , m t , 42 tx 3 6 m 

2£ 8 2£ 9 £ 9 e 10 e 10 

*3(£,t) = ^ + ^-^-^%- (A3) 

£ 7 2 ^7 2 2 £ 8 y J 

4 ' 2£ 2 2£ 2 £ 4 

hh i h a P 2 T 9h 2 P 2 T h\T hh a P 2 h 2 P 2 



2^ 2C 2£ 4 £ 2 R 4 ^ R 4 

U3 TTSS P 2 10 hh^P 2 5u 3tt( P 2 15 h 2 P 2 

+ 2i 2 R 4 e 3 R 4 £ 3 R 4 + C 4 R 4 

15 u 3tt P 2 hfi£ u 3 ttS 2ii 3xx u 3? 



2 



ri? 4 c^ 4 2^ i? 4 cr 4 cr 4 2^ 

15u 3fK -P 4 54m 3 « ^ 4 210^3; P 4 360 u 3 P 4 3u 3;K ^ 2 



*5&t) = 



-I lii I IS _ I 

2^ 3 e e e e 

79m 2{ , P 2 u l££ P 2 427 u 3i P 2 13m t P 2 456 u 3 P 2 

2t 4 ~ + ^ i ~ + ^e e i 6 + 

48 u x P 2 3m 3 ^ 25m 3 ^ u u 34 u 3 2u x 

f 2^ 4 2^ 5 f e 6 e 6 

/i 2 P 2 T 2 M6 ^ 3 P 2 2^ 6 ^ 2 P 2 2m 6t £ 2 u 2rT f 
2£ 4 R 4 R 4 R 4 R 4 

2Au 3 P 2 u 3i Au 3 4ui 

^6 U 8^ + ^-^7T + ^7r + 2u2 (A5) 
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